Zero-bias conductance peak split in d-wave superconductors: 
Derivation of the universal magnetic field dependence 
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The zero-bias conductance peak in d-wave superconductors splits in an applied magnetic field. 
In this work, the experimentally observed universal relation 8 oc y/~Bo for strip-shaped samples is 
derived analytically based on the long-ranged current contributions from Abrikosov vortices inside 
the sample. The result is in full agreement with observed key properties, and features such as 
hysteresis effects are made accessible. Employing a magnetically induced additional order parameter 
is not necessary for the physical explanation of the universal relation. 
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It is well known, that the surface of a material can dis- 
play a different physical behavior than the bulk. Nev- 
ertheless, a thorough analysis of these surface effects 
sometimes reveals important and unique information 
about the bulk system itself. An amazing example are 
fermionic states, commonly referred to as Andreev bound 
states, at the surface of unconventional superconductors 
0, II S S, H 0, 0, 0, EH E3 • Living on a surface sheet of 
a few coherence lengths only, these Andreev bound states 
are responsible for significant zero-bias anomalies in tun- 
neling experiments: They appear in both singlet high-T c 
cuprates [lJ,l2jj3j,J4,|9,|6|]and triplet superconductors such 
as Sr2Ru04 [7|, |8|, [9|, [l0( and they are also expected to 
form in non-centrosymmetric superconductors of mixed 
parity like CePtsSi [ll[. For the fundamental question of 
identifying the superconducting gap function of a given 
bulk material, the surface bound states have proven to be 
of great importance y, y, U, la] • Their mere existence is 
a clear indication, that the pairing interaction underlying 
the superconductivity in the sample is not of the conven- 
tional isotropic s-wave type. Moreover, detailed mea- 
surements of the zero-bias anomalies for different surface 
orientations provide helpful information about the actual 
gap structure in momentum space [H,0] ■ Apart from their 
existence and orientational dependence, another valuable 
source of information is based on the behavior of the zero- 
bias anomalies under the influence of an applied magnetic 
field. As an example, the chirality of a p-wave supercon- 
ducting gap function might be detected due to the unique 
response of the surface bound states, reacting differently 
on magnetic fields of the same magnitude, but antiparal- 
lel orientation For a d- wave superconductor such 

as YBa2Cu307_a;, which is in the focus of this work, the 
conductance peak at zero bias is known to split into two 
peaks at finite energies ±d when a magnetic field is ap- 
plied parallel to the c axis [3, EH, EE E3 • Whereas the 
experimental evidence for the splitting is clear and fully 
established already for quite a long time, the question 
about its physical origin is still open. The main reason 
for this is the existence of different physical mechanisms, 



which are capable of explaining this behavior in general. 
First of all, a natural explanation for the splitting is a 
Doppler shift of the Andreev bound state energies due to 
the superfluid velocity of the condensate at the surface 
[ill. It is important to realize, however, that due to the 
Doppler shift not only Meissner currents screening the 
applied magnetic field contribute to the total splitting, 
but also Abrikosov vortices from inside the superconduc- 
tor, since the currents surrounding them are long-ranged 
1^ . 20j | . A different mechanism, which would also lead to 
a splitting of the zero-bias conductance peak, is the ex- 
istence of a small minority order parameter component 
of the id xy or is type in addition to the main d x 2_ y 2 
order parameter. For this scenario, the magnitude of 
the observed split should be directly proportional to the 
magnitude of the minority order parameter component 

H- 

A very important quantitative observation, which can 
be of great help in identifying the origin of the zero- 
bias conductance peak splitting, has been established 
in a series of experiments done by the Deutscher group 
USUI. Perfor ming tunneling measurements on vari- 
ous narrow, strip-shaped samples of YBa2Cu307_a;, they 
found out, that the zero-bias conductance peak split for 
decreasing magnetic fields obeys the law 



5 = a\J Bq. 



(1) 



Here, Bq denotes the magnitude of the applied magnetic 
field, and a can be determined from experiment. Inter- 
estingly, the asymptotic behavior of all the samples is 
universally described by Eq. with a w l.lmeVT 5 

There 

is a low-field deviation, which seems to 
correlate to the individual doping level of the sample, but 
already for intermediate magnetic fields the above rela- 
tion universally holds true, independent of the specific 
doping. The width of the measured samples varies be- 
tween 600 A and 3200 A. Amazingly, it does not have any 
influence, either. The only exception reported sofar is a 
sample, which is considerably smaller in width, and fulfils 
Eq. (FT]) with a « 0.65meVT' [13 . Up to now, the ex- 



2 



perimental observation of the universal relation has often 
been linked to a theoretical model proposed by Laughlin 
25]. This model assumes the formation of an additional 
id xy order parameter in the bulk of the superconductor, 
which should be induced by an applied magnetic field. In 
this scenario, the magnitude of the induced order param- 
eter, and thus also the corresponding energy split, would 
obey a relation of the same non- linear type as Eq. (fTJ), 
namely S oc \/Bq. It is important to note, however, that 
Laughlin's model does not at all include the contribution 
of Abrikosov vortices in the sample. Since the universal 
relation Eq. ([I]) is measured for applied magnetic fields 
of up to 16 T, which are generally much higher than the 
lower critical field, there are undoubtedly huge numbers 
of Abrikosov vortices inside the samples. Therefore, it is 
essential to properly take into account their contribution 
to the superfluid velocity at the surface. In the following 
it will be clarified, that the presence of Abrikosov vor- 
tices already results in the universal relation Eq. (fTJ) and 
explains its key features. In particular, it is not neces- 
sary to rely on an additional magnetically induced order 
parameter. 

As a starting point for the derivation, we consider a 
d-wave superconducting sample in an applied magnetic 
field Bq 3> B c i, so that Abrikosov vortices are expected 
inside naturally. It is beyond the scope of this work 
to calculate the positions of all these vortices ab initio. 
Rather, it is assumed for simplicity, that the Abrikosov 
vortices form a regular lattice Ao, which is situated at an 
offset do from the boundary. A stability condition ensur- 
ing the physical relevance of this model is employed at a 
later stage. Characteristic samples of such an Abrikosov 
lattice are depicted in Fig. [H a)-d), for two different lat- 
tice types, namely quadratic and triangular ones, at dif- 
ferent orientations. The lattices in Fig. [TJ a) and b), 
have a connection between nearest neighbours, which is 
parallel to the boundary, whereas those in c) and d) are 
rotated. It is important to note, that we may find lines 
parallel to the boundary (marked as dashed) , which inter- 
sect all vortices in such a way, that two adjacent vortices 
on such a line always have the same distance s. More- 
over, the spacing between these lines is equidistant and 
denoted by d in the following. Obviously, the area of 
a unit cell of the lattice is given by the product s ■ d, 
whereas the ratio q\ = d/s depends on both the lattice 
type and its orientation, which will have consequences on 
the resulting zero-bias conductance peak split as shown 
below. 

Abrikosov vortices are surrounded by long-ranged su- 
percurrents, which decay on the lengthscalc of the mag- 
netic penetration depth A. For the purpose of this work, 
these supercurrents are well approximated by the vortex 
solution of London electrodynamics, since the exact be- 
haviour in the core region is of no relevance here. In the 
following, the real gauge is employed, so that a phase 
gradient has already been absorbed into the magnetic 
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FIG. 1: (color online) Sketch of different regular vortex lat- 
tices situated at an offset do from the boundary of the su- 
perconductor. Each of the lattices can be characterized by 
specific lines, which are parallel to the boundary (dashed). 
Along such a line, adjacent vortices have a distance s, and 
two neighbouring lines themselves are separated by the dis- 
tance d. In units of the intervortex distance between nearest 
neighbours, we have a) quadratic lattice: s = d = 1, b) trian- 
gular lattice: s = 1, d = s/3/2, c) rotated quadratic lattice: 
s = \/~2, d = 1/V2, d) rotated triangular lattice: s — v3, 
d = 1/2. The ratio q\ — d/s depends on both lattice type 
and orientation. 



vector potential A, which itself is then proportional to 
the superfluid velocity. The London solution of a single 
free Abrikosov vortex is given by 

a 1 s *o „ f \r~r v \ \ (x - x v )y - (y - y v )Z 

(2) 

where $0 = hc/2e is the flux quantum, K\ a modified 
Bessel function of first order, and r = (x, y) the space 
coordinate. Accordingly, r v = (x v ,y v ) denotes the posi- 
tion of the vortex. The correct vortex solution for a ge- 
ometrically restricted area has to fulfil Neumann bound- 
ary conditions, i.e., the superfluid velocity and thus the 
magnetic vector potential must be parallel to the bound- 
aries. Generally, the free vortex according to Eq. ([2]), 
which represents a circular vector field around r„, does 
not accomplish this condition. For some simple geome- 
tries, however, it is possible to construct the solution by 
placing virtual vortices and antivortices outside the su- 
perconducting area, which is analogous to the technique 
of mirror charges commonly known from classical elec- 
trostatics. For the concrete case of the vortex lattice Ao 
in a superconducting strip of width L considered here, 
the correct solution can be generated by adding infinitely 
many virtual lattices A„^o> consisting of (anti) vortices 
for n even (odd). A sketch of the situation is shown in 
Fig. [21 Obviously, this arrangement of vortices has the 
required symmetry of exchanging vortices and antivor- 
tices when mirrored at the boundaries x = or x = L, 
respectively, resulting in the correct implementation of 
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FIG. 2: (color online) Sketch of a regular vortex lattice Ao 
inside the d-wave superconducting strip (blue) of width L. 
The orientation of the d wave is indicated by the lobes. To 
ensure the correct implementation of electrodynamic bound- 
ary conditions, an infinite number of virtual mirror lattices 
A n ^o has to be taken into account. For odd n, A n consists of 
antivortices (red). 



the Neumann boundary conditions. Thus, the magnetic 
vector potential due to the Abrikosov vortices can be 
written as 



E 



E 
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(3) 



Note, that this contribution from the Abrikosov vortices 
corresponds to a magnetic field Ba = V x Aa that van- 
ishes along the boundaries per construction. An external 
magnetic field Bo applied along the z-axis has to be taken 
into account separately by the well-known Meissner so- 
lution of London theory: 



cosh(z/A-£/2A) 
BM(r) - B ° cosh(L/2A) Z 
sinh(a;/A - L/2X) , 

Aoo , , r — y • 



A M (r) 



cosh(L/2A) 



(4a) 
(4b) 



Now that both the contributions of Abrikosov vortices 
and Meissner currents are found, the total magnetic vec- 
tor potential is just their superposition: A = Am + Aa- 
In the following, the zero-bias conductance peak split 
5 = |«f|A|, which is expected at the boundary due to 
the Doppler shift of the superfluid velocity, will be calcu- 
lated. Note, that the relevant quantity for the tunneling 
measurements is S, which is spatially averaged along the 
boundary. The Meissner contribution is constant along 
the boundary anyway, and the result at x — directly is 



A M = -A£ tanh(L/2A)y. 



(5) 



The contribution of the vortex lattice Ao and its mirror 
lattices A„^o is periodic on the length s (cf. Fig. []}. 
Therefore, their spatial average along the boundary is 



given as 



Aa = - 



dyA A (x = 0,y). 



(6) 



It it possible to integrate out the y-dependence analyti- 
cally, which leads to 
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(7) 



n=Q i„£A„ 



The exponential term is the contribution of all Abrikosov 
vortices with the same a^-coordinate, corresponding to a 
whole dashed line in Fig. [TJ Indexing these positions x v 
according to x v — nL + md + do with m = 0, 1, m c — 
(L — 2do)/d leads to 
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(-l) n e~ nL/x e- md/x y. (8) 



Analytic evaluation of both geometric series results in 
- $ sinh(L/2A + rf/2A-rf /A) „ 



2s cosh(L/2A)sinh(d/2A) 



(9) 



At this stage, it is essential to employ an additional con- 
dition, which ensures the stability of the vortex lattice. 
For a specific vortex of the lattice to remain stable, the 
current density generated by all other vortices, antivor- 
tices and the Meissner screening current should vanish at 
the corresponding vortex position. After a rather lengthy 
but straightforward analysis similar to the one above, the 
following condition can be derived 



£o = . sinh(d/ 2A) 
2s °cosh(d/2A 



do/A) : 



(10) 



which guarantees the stability of all vortices simultane- 
ously, i.e., the whole vortex lattice is stable. Plugging 
the stability condition Eq. (fTOj) into Eq. @ yields 



A A = XB [tanh(L/2A) + tanh(d/2A - d /X)] y. (11) 

Note, that the first term, which depends on the width 
L and is linear in the applied magnetic field amplitude 
Bo, exactly cancels the Meissner contribution Am of Eq. 
([5|). Only the remaining second term of the vortex lat- 
tice contributes to the total superfluid velocity, so that 
eventually 



'-v F XB |tanh(d/2A - d /X)\ 



(12) 



It is important to realize, that this remaining term is not 
linear in Bo, since the parameters d and do characterizing 
the vortex lattice scale as a function of the magnetic field 
themselves. However, because d/X,do/X <C 1 already for 
intermediate fields Bo, Eq. (jlpp may be expanded in 
these small parameters, yielding d = (gA^o/^o) 5 - Ac- 
cordingly, the final result for the zero-bias conductance 
peak split due to the Abrikosov vortices is 



2JW, 



--Vf/aVBo 



(13) 
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FIG. 3: (color online) Zero-bias conductance peak split 8 as 
a function of the applied magnetic field Bo. Main figure: Ex- 
perimental data for sample S4 of Ref. [2J] in decreasing fields 
(red points) and theoretical fit according to Eq. (113[) with 
/a = 0.0845, which follows e.g. for a triangular Abrikosov vor- 
tex lattice with an offset of do/d = 0.5909. Inset: Experimen- 
tal data for the 600 A sample of Ref. Fig. 3, in increasing 
and decreasing fields (blue points). Curves are sample fits for 
a triangular vortex lattice at an offset of do/d = 0.6554 (in- 
creasing field) and do/d = 0.5923 (decreasing field). 



This work has shown, that the long-ranged current 
contributions from Abrikosov vortices inside a e?-wave 
superconductor are directly responsible for the univer- 
sal magnetic field dependence of the zero-bias conduc- 
tance peak split observed experimentally. Based on this 
derivation, key properties are fully explained and other 
observed features are accessible qualitatively. Deviations 
of the average Abrikosov vortex density from that of a 
regular lattice, e.g. due to strong pinning, expected to 
result in deviations from Eq. |T]), are beyond the scope of 
this work. Nevertheless, the zero-bias conductance peak 
split for arbitrary vortex distributions may be derived 
similarly as presented here in a straightforward way. 

I acknowledge valuable discussions with N. Schopohl, 
T. Dahm, and M. Sigrist, and I am grateful to R. Beck for 
providing me with the experimental results. This work 
was financially supported by the Swiss Nationalfonds and 
the NCCR MaNEP. 
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with the dimensionless number f\ = ^JqX 

The derived Eq. (|13[) is in complete agreement with 
the experimentally established universal relation Eq. fl} . 
Besides natural constants, the prefactor consists of the 
Fermi velocity v f as a material parameter and the factor 
/a, which is determined by the Abrikosov vortex lattice 
type and offset. Reasonable assumptions for the vortex 
lattice geometry allow a perfect fit to experimental data 
by Beck et al. [H, Hi], cf. Fig. H Apart from the 



quantitative explanation of the universal magnetic field 
dependence, the Abrikosov vortices might also provide 
some qualitative insight to the commonly observed hys- 
teresis effects of the splitting. The number of Abrikosov 
vortices in a given sample is higher, if the magnetic field 
Bo is accessed downwards from higher magnitudes com- 
pared to the case, where the magnetic field has been in- 
creased. For equal shape of the vortex lattice and equal 
intervortex distance, this directly yields a larger offset 
parameter do/d at increasing magnetic fields and corre- 
spondingly a larger splitting, cf. the inset of Fig. [3] 
Finally, it should be pointed out that the resulting Eq. 
(jT5]l does not depend on the width L of the superconduct- 
ing strip any more, as experimentally observed. Albeit 
speculative, even the one exceptional sample of differ- 
ent a and much smaller width might be explained. It 
is reasonable to assume, that a phase transition of the 
vortex lattice structure is induced geometrically by the 
further narrowing of the strip width, here. For example, 
an induced phase transition from triangular to rotated 
triangular shape (cf. Fig. [T|) would lead to a different 
value of q\ and result in an amplitude a different by a 



factor of w 0.58, which is close to the reported value 

of ^ 
Ui i.i 



0.59. 
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